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Abstract: Simulations of geosynchrotron radio emission from extensive air showers performed with the 
Monte Carlo code REAS1 used analytical parameterisations to describe the spatial, temporal, energy 
and angular particle distributions in air showers. The successor REAS2 replaces these parameterisations 
with precise, multi-dimensional histograms derived from per-shower CORSIKA simulations. REAS2 
allows an independent selection between parameterisation and histogram for each of the relevant particle 
distributions, enabling us to study the changes arising from using a more realistic air shower model in 
detail. We describe the new simulation strategy and illustrate the effects introduced by the improved air 
shower model. 
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Introduction 

From results of the LOPES experiment [1], radio 
emission from cosmic ray air showers is known to 
be dominated by a geomagnetic emission mech- 
anism that can be described with the "geosyn- 
chrotron model" [2]. In the geosynchrotron pro- 
cess, relativistic secondary shower electrons and 
positrons are deflected in the earth's magnetic 
field, thereby giving rise to strongly pulsed, co- 
herent radio emission in the frequency range from 
~ 10 to 100 MHz. In the recent years, the 
geosynchrotron model had evolved from analytic 
frequency-domain calculations [4] to time-domain 
Monte Carlo simulations based on an analytical de- 
scription of the underlying extensive air shower, 
as implemented in the REAS1 simulation code 
[5]. Many general properties of the radio emis- 
sion have been predicted since using REAS1 [6]. 
The successor to this code, REAS2, now features 
an implementation of the geosynchrotron model 
no longer based on analytically parameterised air 
shower properties, but using realistic, per-shower 
CORSIKA [3] simulations to model the spatial, 
temporal, energy and angular distributions of the 
shower particles. In this article, we describe the 
changes arising from the much more realistic air 
shower model and illustrate the potential of the 



new simulation code for advanced studies of the 
radio pulse shape. 

REAS2 air shower model 

While the electromagnetic emission model has not 
changed between REAS1 and REAS2, REAS2 
simulations are now based on very detailed parti- 
cle information derived with CORSIKA on a per- 
shower basis. For each individual air shower, 
CORSIKA writes out separate information for 
electrons and positrons sampled in (usually) 50 
layers between the point of first interaction and the 
observer position. Each layer encompasses 

• one three-dimensional histogram of 

1. particle arrival time relative to that of 
an imaginary primary particle propa- 
gating with the speed of light from the 
point of first interaction 

2. lateral distance of the particle from the 
shower core 

3. particle energy 

• and one three-dimensional histogram of 

1 . angle of the particle momentum to the 
shower axis 



REAS2 GEOSYNCHROTRON SIMULATIONS 



2. angle of the particle momentum to the 
(radial) outward direction 

3. particle energy. 

These histograms give REAS2 access to a true 
four-dimensional distribution of particles in atmo- 
spheric depth, arrival time, lateral distance and en- 
ergy, and in addition describe the angular distri- 
bution of particle momenta as a function of par- 
ticle energy and atmospheric depth. The chosen 
separation of the distributions into two histograms 
ensures that the necessary amounts of data can be 
handled on standard PCs while making approxima- 
tions of only minor significance for the simulation 
of the radio signal. (The most important drawback 
of this scheme is the loss of information on az- 
imuthal asymmetries in the air shower.) Naturally, 
effects associated with air showers induced by dif- 
ferent types of primary particles can be analysed in 
detail with this simulation strategy. The longitudi- 
nal evolution of the air shower is sampled on an ad- 
ditional, finer grid (usually) spaced with 5 g cm -2 
distance. 

REAS2 vs. REAS1 results 

Incorporating the new CORSIKA-based air shower 
model and some additional enhancements, the 
REAS2 code provides a much more realistic deter- 
mination of geosynchrotron radio emission than its 
predecessor REAS 1 . A particular merit of the cho- 
sen approach is that the transition from the REAS 1 
to the REAS2 air shower model can be performed 
in a gradual fashion, switching the different parti- 
cle distributions (spatial, temporal, energy, and an- 
gular) from parameterised to histogrammed one at 
a time and analysing the changes arising in the ra- 
dio signal. The corresponding effects and techni- 
cal details have been discussed elsewhere [8]; here 
we only compare the end result of REAS2 simu- 
lations with those of REAS1 simulations for the 
typical reference case of a vertical 10 17 eV proton- 
induced air shower. As presented in Fig. 1, the ra- 
dio pulses calculated with REAS2 for an observer 
to the north of the shower core show only moderate 
changes in comparison with the REAS 1 -generated 
pulses. In particular, the pulses close to the core 
(75 m corresponds to the typical lateral distance 
in the LOPES experiment) get considerably nar- 



rower, caused by the narrower arrival time dis- 
tributions provided by CORSIKA in comparison 
with the parameterisation used in REAS1. Con- 
sequently, the frequency spectrum of the emission 
close to the core gets much flatter for the REAS2- 
calculated pulses (Fig. 2). Further away from the 
core (525 m corresponds to the distance range of 
interest for larger scale radio antenna arrays), the 
amplitude drops by a factor of ~ 2, mainly as 
a consequence of the much broader angular dis- 
tribution of particle momenta derived from COR- 
SIKA. Interestingly, the overall field strength in 
the frequency band used by the LOPES experiment 
(40 to 80 MHz) does not change significantly be- 
tween REAS1 and REAS2 (Fig. 2). A more sig- 
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Figure 3: Comparison of REAS1 and REAS2- 
simulated radio pulses at 75 m and 525 m west of 
the shower core. The bipolar pulses present in the 
REAS 1 simulations are no longer present. 

nificant change becomes visible when studying the 
changes to the radio pulses for an observer west of 
the shower core, as depicted in Fig. 3. First, a sig- 
nificant drop in pulse amplitude close to the shower 
core can be identified. This leads to a pronounced 
east-west versus north-south asymmetry in the ra- 
dio "footprint", even for vertical air showers (cf. 
Fig. 4). At larger distances, a qualitative change 
in the pulse shape takes place: while REAS1- 
generated pulses in this region showed bipolar 
structures, the REAS2-calculated signals become 
universally unipolar. The bipolar pulse shapes in 
the REAS1 calculations can be considered arti- 
facts of over-simplified particle distributions (e.g., 
in the momentum angles), and the REAS2 results 
describe the geosynchrotron emission much more 
realistically. Another important characteristic of 
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Figure 1: Comparison of REAS1 and REAS2- 
simulated radio pulses at 75 m and 525 m north 
of the shower core. 
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Figure 2: Comparison of frequency spectra for the 
REAS1 (thin) and REAS2 (thick) simulated radio 
pulses. 
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Figure 4: Contour plots of REAS1 (upper) and REAS2 (lower) air shower emission at v = 60 MHz. The 
columns (from left to right) show the total field strength, the north-south polarisation component and the 
east-west polarisation component. The vertical polarisation component (not shown here) does not contain 
any significant flux. Contour levels are 0.25 /iV vaT 1 MHz -1 apart in E u , outermost contour corresponds 
to 0.25 /iV m _1 MHz -1 . White centre region has not been calculated. 



geosynchrotron radiation is its mostly linear po- 
larisation. The comparison of the individual po- 
larisation components depicted in Fig. 4 demon- 
strates that the polarisation characteristics are iden- 
tical between REAS1 and REAS2. At the same 



time, the contour plots illustrate once more the 
newly arising east-west versus north-south asym- 
metry and confirm that the absolute field strengths 
in the centre of the LOPES band (60 MHz) do not 
change considerably. 
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Figure 5: Contribution of different longitudinal 
shower evolution stages to the radio pulse at 75 m 
north from the shower core. 

Pulse shape analyses 

The highly detailed CORSIKA-based air shower 
model implemented in REAS2 allows advanced 
studies of the radio pulse shape. A particularly in- 
teresting question is how different phases in the 
longitudinal development of the air shower con- 
tribute to the radio pulses, as illustrated in Figs. 5 
and 6. Close to the shower core, signals from all 
over the longitudinal shower evolution arrive ap- 
proximately simultaneously at the observer. The 
pulse shape close to the shower core thus gives a 
direct estimate of the overall particle arrival time 
distribution. (This could change once the refrac- 
tive index profile of the atmosphere is taken into 
account.) At larger distances, geometrical time de- 
lays become important; the pulse shape thus pro- 
vides direct information on the shower evolution 
profile. Another interesting result is that the emis- 
sion is dominated by the shower maximum (here 
at 640 g cm~ 2 ) and the stage shortly before. The 
information content of the radio pulses can be ex- 
ploited to estimate the primary particle energy and 
type from radio measurements on a shower-to- 
shower basis [7]. Analyses how different particle 
energy ranges or radial distance ranges contribute 
to the radio signal have also been performed [8]. 

Conclusions 

With REAS2, a sophisticated Monte Carlo im- 
plementation of the geosynchrotron model based 
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Figure 6: Contribution of different longitudinal 
shower evolution stages to the radio pulse at 525 m 
north from the shower core. 

on a realistic, CORSIKA-based air shower model 
is now available. The transition from REAS1 
to REAS2 has been carried out in a controlled 
way, and the changes arising are well-understood. 
In spite of the major model improvements, the 
changes are only moderate, in particular in the fre- 
quency range of current experiments. REAS2 can 
be used for in-depth studies of the information con- 
tent of geosynchrotron radio pulses and as such is 
a powerful tool to unlock the full potential of the 
radio technique for cosmic ray measurements. 
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